Carrier induced ferromagnetism in diluted local-moment systems 
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The electronic and magnetic properties of concentrated and diluted ferromagnetic semiconductors 
are investigated by using the Kondo lattice model, which describes an interband exchange coupling 
between itinerant conduction electrons and localized magnetic moments. In our calculations, the 
electronic problem and the local magnetic problem are solved separately. For the electronic part an 
interpolating self-energy approach together with a coherent potential approximation (CPA) treat- 
ment of a dynamical alloy analogy is used to calculate temperature-dependent quasiparticle densities 
of states and the electronic self-energy of the diluted local-moment system. For constructing the 
magnetic phase diagram we use a modified RKKY theory by mapping the interband exchange to an 
effective Heisenberg model. The exchange integrals appear as functionals of the diluted electronic 
self-energy being therefore temperature- and carrier-concentration-dependent and covering RKKY 
as well as double exchange behavior. The disorder of the localized moments in the effective Heisen- 
berg model is solved by a generalized locator CPA approach. The main results are: 1) extremely low 
carrier concentrations are sufficient to induce ferromagnetism; 2) the Curie temperature exhibits a 
strikingly non-monotonic behavior as a function of carrier concentration with a distinct maximum; 
3) Tc curves break down at critical n/x due to antiferromagnetic correlations and 4) the dilution 
always lowers Tc but broadens the ferromagnetic region with respect to carrier concentration. 

PACS numbers: 75.50.Pp, 71.23.-k, 75.30.Hx, 85.75.-d 
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I. INTRODUCTION 



Diluted Magnetic Semiconductors (DMS) based on III- 
V and II- VI semiconductors doped with magnetic impu- 
rities such as Mn have attracted great experimental and 
theoretical interest in recent years due to their poten- 
tial in spintronics applications 1 . One of the most widely 
studied DMS is Gai-^Mn^As of which Curie tempera- 
ture as high as Tc = 173K has been observed!. These 
materials indeed exhibit a very striking correlation be- 
tween the transport and magnetic properties. Ideally, 
each Mn dopant atom represents an acceptor that in- 
troduces a local spin and a hole carrier. The ferromag- 
netism is driven by a charge-carrier mediated mechanism 
as a consequence of an interband exchange interaction 
between the localized magnetic moments and the car- 
rier (hole) spins. Empirically, one finds that the hole 
density in DMSs is significantly smaller than the dop- 
ing concentration due to carrier compensation etc. In 
addition, the measured saturated magnetization is ob- 
served to be much less than the nominal concentration 
of Mn ions 2 . Experimentally, the remanence magnetiza- 
tion and the Curie temperature Tc are both enhanced 
by post-growth annealing of the samples, which changes 
positions of defects and the hole concentration 3 -^. So 
the magnetism in DMSs is heavily affected by material 
issues, e.g. very disordered impurity-lattice-location, car- 
rier compensation and the coupling mechanism between 
the localized Mn spins and the itinerant holes, as well as 
a complicated band structure. Therefore, it is impossible 
to study these materials without making some approxi- 
mations. 

There are various theoretical approaches attempting 



to understand DMS physics. One class of models ne- 
glects interband coupling effects and only investigates 
the disorder of the localized moments^&i^. Some mod- 
els introduce the interband coupling between the itin- 
erant carriers and the localized moments and study 
the effect of the charge carriers on the magnetic phase 
diagran^^ii^iidSiL 17 .. If the ferromagnetic or- 
dering is mediated by 'free' carriers, its origin is as- 
sumed to be of RKKY type 18 . However, further re- 
sults show that RKKY interaction between the local- 
ized spins is insufficient 3 - ^ 19 ! 20 ! 21 ' 22 . Korringa-Kohn- 
Rostock CPA (KKR-CPA) calculations^! show that 
the ground state of Gai-^Mn^As is half-metallic. The 
k • p approac h 25 ' 26 i 27 ' 28 uses the p — d exchange inter- 
action between Mn and hole spins, in which the states 
near the Fermi energy are mainly assumed to have 
the character of the host semiconductor valence band. 
Moreover, dynamical mean-field theory (DMFT) 29,30,31 
is used to investigate the carrier spin polarization 
of the impurity band and the Curie temperature of 
Gai-a-Mn^As. In addition to these theories, ab-initio 
approache a 32 ' 33 '^' 3 ^, spin- wave theorie a 36 ' 37 '^' 3 ^, Monte 
Carlo (MC) simulations^^ -^-^ 3 - have been applied to 
understand the physics of DMSs. 

However, both interband-coupling effects and the dis- 
order effects on the ferromagnetic phase transition have 
not been completely understood so far. Instead, the dif- 
ferent theoretical models have tended to concentrate on 
different aspects of the problem. In this paper, we will 
use the Kondo lattice model (KLM) to investigate effects 
of the interband coupling and the disorder on the mag- 
netic and electronic properties of diluted local-moment 
systems such as DMS. The positional disorder of the mag- 
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netic ions and the interband exchange coupling of charge 
carriers to the disordered local moments are treated on 
the same level and in the same theoretical framework. 
Here, we neglect the Coulomb disorder potential due to 
the substitutional Mn impurities. As a model study, the 
concentration of magnetic ions is not limited to x ~ 0.05 
as in real DMS. We investigate instead the whole concen- 
tration range x = 0- • • 1. The paper is expected to pro- 
vide the qualitative trends of the interconnected effects of 
interband coupling and disorder on the magnetic proper- 
ties of diluted ferromagnetic semiconductors. Quantita- 
tive agreement with experiment will require treating both 
low magnetic ion concentration and defect correlations in 
a realistic way 2 ". 

The article is organized as follows. The theoretical 
methods are described in section [TTJ Section IIIII is con- 
cerned with the numerical results and a discussion of 
them. Section IIVI concludes the article with a summary. 

II. MODEL 

The Kondo lattice model (KLM) is surely one of the 
most intensively discussed models over the past years in 
the field of magnetism. It describes the interaction be- 
tween itinerant carriers (spin cr^) in a partially filled en- 
ergy band and localized magnetic moments (spin S^) at 
certain lattice sites: 

H s f — — °» " Si , (1) 

i 

where the index i refers to the lattice site and J is the 
coupling between the itinerant carriers and the localized 
moments. The respective model Hamiltonian consists of 
two parts: 

H = H S + H sf . (2) 

H s describes itinerant carriers in an energy band of finite 
width 

ija 

where is the hopping integral, cj ff (c !ff ) creates (anni- 
hilates) a carrier on the site i and a = \, j is the spin pro- 
jection. According to the sign of the exchange coupling 
J, a parallel (J > 0) or an antiparallel (J < 0) alignment 
of itinerant and localized spin is favored with remarkable 
differences in the physical properties. We restrict our 
considerations in the following to the J > case, some- 
times referred to as 'ferromagnetic Kondo-lattice model 
(FKLM)' (s — /, s — d model) or, in the strong-coupling 
regime, as the 'double-exchange model'. The so-called 
p — d exchange, not considered here, would require J < 
and different, significantly greater, absoluted values^. 
At first glance, one does not expect dramatic changes 
of the collective magnetic order since in the first order 



the features follow J 2 . In the following, the charge carri- 
ers will be itinerant electrons. The generalization to the 
case of holes is straightforward. 

The KLM refers to magnetic materials that get their 
magnetic properties from a system of localized magnetic 
moments being indirectly coupled via interband exchange 
to itinerant carriers. The magnetic properties of the 
FKLM, therefore, strongly correlate with the electronic 
properties of the itinerant electron subsystem. Thus, the 
evaluation of the FKLM includes the solution of two in- 
terconnected partial problems, mainly an electronic prob- 
lem and a magnetic local-moment problem. In addition, 
for the diluted systems, the inherent positional disorder 
of the magnetic ions does not only break the translational 
symmetry of the crystal lattices but also influences the 
electrons distribution. The disorder of the magnetic-ion 
position and the interband coupling greatly complicates 
the theoretical description of the FKLM, too. Therefore, 
the evalution of the FKLM has to deal with the conduc- 
tion band disorder and on-site local-moment disorder. So 
it is really difficult to investigate all those factors instan- 
taneously as the same level. In this paper, accordingly we 
apply a strategy of self-consistent calculation - the pro- 
cedure is mainly divided into two parts: I) the electronic 
part in which the local-moment magnetization is fixed as 
parameter; 2) the local magnetic part in which the solu- 
tion of the electronic part leads to an effective exchange 
integral for the interaction between the local spins. The 
following sections III Al III Bl III CI and III PI show the theo- 
retical methods and the section Hi El describes the whole 
loop of self-consistent calculation. In the last analysis, 
the model description refers to the mechanism of carrier- 
mediated ferromagnetism as it is the case in diluted III- 
Mn-V and n-doped II-Mn-VI samples. 

A. The concentrated Kondo-lattice model 

The model Hamiltonian © provokes a rather sophis- 
ticated many-body problem. Using second quantization, 
the interaction term |T]) can be rewritten as: 

where z a = 5^ - 5 al , SJ = Sf + iz a S^ and n i<7 = c\ a c la 
is the occupation number operator at site i. The first 
term on the r.h.s of (H|) shows an Ising-like interaction 
between the z-components of the spin operators and the 
second term describes spin-exchange processes. 

The electronic part of the many-body problem of the 
above Hamiltonian is solved as soon as one gets the re- 
tarded single-particle Green function G^ a (E): 

Gf 3a {E) = {{c l<y ;c] a )) E = ±J2 G t(E)e^ , (5) 

k 
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where the superscript "el" refers to "electron". The Bloch 
energy e(k) is the Fourier transform of the hopping inte- 
gral T lf e(k) = Y~ 4 £ H !] ,< k " and R 



R , R , 



The electronic self-energy Mic a (E) becomes the central 
quantity of the many-body problem of the electronic part. 

For finite temperatures and arbitrary band occupa- 
tions, an exact expression of M-^ a {E) is not available and 
one needs to apply an approximatio n 45 ' 46 . The method 
that we choose is the Interpolating Self- energy Approach 
(ISA)il which results in a wave- vector independent self- 
energy 



M (E) - J -z (S*)\ J2 a ^i,E-jz^)) 



(7) 



The parameters a ai b a are fixed by rigorous high-energy 
expansions to fulfill the first four spectral moments: 

a a = S{S+l)-z a {S z )( Zl7 (S z ) + l) , b a = b^='^ (8) 

and Gq(E) is the free propagator: 



N^E- e(k) 



(9) 



The ansatz ([7]) is chosen to be exact for a maximum num- 
ber of special cases in the low-density limit. Therefore, 
it should represent a reasonable starting point for the 
description of ferromagnetic semiconductors. The ISA 
approach has already been applied in model studies^ and 
to real substances^. For a detailed account of the de- 
coupling procedure and an extensive discussion of the re- 
liability of this self-energy we refer the reader to Ref. 47 . 
M a (E) is the electronic self-energy for the concentrated 
KLM and, in the next subsection, will be used to sim- 
ulate the disorder of the magnetic moments in diluted 
ferromagnetic semiconductors. 



then be derived from the propagator 



R a (E)= Jdu; 



E-LO- Scr(-E) 



(12) 



where is now the electronic self-energy in the di- 

luted system and po(x) the Bloch-density of states of 
the non-interacting carriers. For the determination of 
the decisive self-energy, one can use a standard CPA 
formalism 47 : 







ca 



Cb 



-X a (E) 



1-# CT (£)(-E CT (£)) 

M a {E)-Jl a {E) 
1 — R<r(E)(M a (E) — Jl a (E)) 



(13) 



where we have chosen £o = 0. 

The configurational averaging, inherent in CPA, 
provides translational symmetry and therewith site- 
independent average spin-dependent occupation num- 
bers: 

+ oo +oo 
— oo — oo 

where f-{E) is the Fermi function, /i the chemical po- 
tential and Pa-(E) the quasiparticle density of states of 
the interacting particle system: 



Pa(E) = --lmR a (E) 



(15) 



Thus, for given magnetization (S z ) and occupation num- 
ber of carrier n = (rtj) + (n^), the many-body problem in 
the diluted system can be solved by the Eqs. ©, (fl~2j) . 
(fT3"]) and (TTJJ in a self-consistent manner. The single- 
particle Green function can be expressed as 



B. Dynamical alloy analogy 



For diluted ferromagnetic semiconductors A CA B CB of 
nonmagnetic ions A' and magnetic ions 'B', where ca 
and cb refer to the concentrations, one should investigate 
the effect of disorder on the electronic self-energy. If the 
atomic level of A sites is 



£Ao — £o , 



(10) 



N^> E-e(k)-E a (E) 



(16) 



The Gf^" (E) describes the electronic properties of the 
itinerant electron subsystem in the diluted ferromagnetic 
semiconductors. In the next subsection, we discuss a 
modified RKKY (MRKKY) theory 4 ^ that will connect 
the magnetic properties and the electronic properties. 



using dynamical alloy analogy (DA A), the atomic level 
of B sites reads as a dynamic atomic energy level: 



e B <r = e +M a (E) 



(11) 



since the local interband exchange H s f acts on the charge 
carriers. Neglecting a Coulomb disorder potential which 
might be important in some circumstances^, (e.g. metal- 
insulator transition), the single-particle properties can 



C. Modified RKKY 

The basic idea of the MRKKY 50 is to map the ex- 
change interaction |T]) on an effective spin Hamiltonian 
of the Heisenberg type: 



H 



(Hsf) 



ij 



(17) 
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by averaging out the conduction electron degrees of free- 
dom ((• • • )( c -*). If one writes H s f as: 

H >t = - -4 E E e ~ iqR< ( s « ■ & )™> ( 18 ) 



AT 

with Pauli spin matrices <r of the carrier-spin operator, 
the averaging (H s f)^ means to retain operator character 
in the subspace of the local spins and to calculate the 
expectation value: 



(4 



k+qer C k£r 



°k+qo- L ko 



(19) 



W has exactly the same structure as the KLM Hamil- 
tonian H Eq. {2|), except for the fact that the local spin 
operators are to be considered as c-numbers for the aver- 
aging (• • • )( c \ therefore not affecting the trace. 3' is the 
corresponding grand partition function. The expectation 
value Eq. (fT9|) does not necessarily vanish for q ^ and 
for a ^ a', as it would do when averaging in the full 
Hilbert space of the KLM. To calculate Eq. (TiT)|) , one 
introduces a proper restricted Green function 



r k,k 



- q (£) = «c kCT ';eL +qff » 



,( c ) 



(20) 



which has the normal definition of a retarded Green func- 
tion, only the averages have to be done in the Hilbert 
space of H' (the details of the averaging procedure being 
in the Ref . The solution of the equation of motion for 
the restricted Green function can be iterated up to any 
desired accuracy, however, at the expensive of higher and 
higher products of local spin operators. To retain the op- 
erator form (fT7|) one has to truncate the infinite series in 
a proper wa y 46 ? 50 . Then, the effective exchange integrals 
of the Heisenberg-Hamiltonian Hf can be written as: 



J?! f = - 



j 2 



47T^ 



Im 



dEf_{E)Gt£{E)Gf£°{E). (21) 



Since ■/?• is a functional of the conduction electrons self- 
energy Y, a (E), the effective exchange integrals receive a 
distinct temperature and band-occupation dependence. 
Neglecting the self-energy T, a (E), J^- leads to the con- 
ventional RKKY interaction. Once getting the exchange 
integrals J?f , one can solve the magnetic properties of 
the KLM from the effective Heisenberg Hamiltonian (TTT)) . 



D. The disordered Heisenberg model 



In this section, a CPA-like approaches! is used to 
study a disordered Heisenberg model with long-range 
exchange. In the following, we simply introduce this 
method and apply it to the diluted spin systems. 

To solve the Heisenberg Hamiltonian 



(22) 



one introduces the retarded magnon Green's function 



Gij(E) - ((Si ; Sj 



\\ ret 
// E 



(23) 



where — Sf ± iSf . For simplicity we suppress the 
index 'eff ' of the exchange integrals. It should be men- 
tioned here that the index i and j in the Eq. (|23| refer 
to all sites of systems and is different from the Hf in the 
Eq. (32). 

For a binary spin systems A CA B CB , the Green's func- 
tions' equation of motion based on the well-known Tyab- 
likov approximation reads^ 



Ej 



if one introduces the 2x2 matrices 



Gf A Gf 



UiGijXj yiGijyj 



(24) 



(25) 



jAA j 



. IB 



and 



Si = 



Xi9i 




Vi9l 



' OJ — UJA 









where Xi + yi — 1 (xj = 1 or if i G A or B, respectively), 
J* Y the exchange integrals between X sites and Y sites 
(X,Y e A or B), a x = (Sf) /a Q (A = A or B, a = 
ca (S A ) + c B (Sf,)), u> 
loa or U)B 



E/2a and the 'local potential' 



m m 
UJ B = C A a A (^J BA ) + CbO-bC^JI, 



J BB ) 

%va I 



(26a) 
(26b) 



Here, the influence of neighbor spins by exchange inter- 
action with site i is taken as an 'effective static field'. 
The random A-type or B-type occupation of any site de- 
termines the 'local potential' loa or lob, respectively. 

Thus the matrix propagator satisfies an equation 
for a problem with only 'local disorder', the matrix JJy 
being independent of the random numbers. Let us now 
extend the theory to the special case of a diluted mag- 
netic systems we are exclusively interested in here. For 
this purpose, we assume that A is a nonmagnetic ion, 



J 



BA 



0, 



i.e., the exchange integrals J AA = J AB 
the averaging (S\) = and (S B ) — (S z ). In addition, 
J BB describes the exchange integral of magnetic ion B 

and therefore J BB = J^j^ that can be calculated by Eq. 

{H}. Thcn . on e gets gf = 0, G AA = and 



eff 

ij 

eff 



9? f (k 
Vi9? = 



~ E ^im^mj) ! ( 27 ) 
m 

Viiu'-CB^J&r 1 . (28) 
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Here G e J f = Off and J = E/{2(S Z )). Thus, the com- 
plete averages lead to (G) = (G e ^). 

Following the locator approach in ReL 7 -, one introduces 
the renormalized locator: 

l=(9? f /(l-9? f U )), (29) 

where Uq is the diagonal element of the renormalized in- 
teractor, which describes the interaction in the effective 
medium 7 . Thus the CPA self-consistent condition is ex- 
pressed as£2 



eff 



(30) 



where G^. = [7 1 



U Q + JeffQi.)}- 1 and J e// (k) is the 
Fourier transform of the exchange interaction . The 
renormalized locator can be written as 



7 (u/) = 



CB 



eff 



(31) 



Once one has found the renormalized locator 7 from 
self-consistent CPA equations Eq. (|3"0|) and Eq. (|3"Tj) , the 
magnetization can be expressed by the Callen equation 53 : 

czx _ (s-m + $) 2S+1 + (i + s + j>)<i> 2s+i 

^ ; (1 + $)2S+1 _$2S+1 ^ 

where the average magnon number can be calculated by 

D(u/) 



du;'- 



(33) 



e 2{S*)w'/k B T _ I ' 

Imry(uj') / (t:cb) is the magnon spectral 



where D(u>') 
function. 

Now, for a given temperature and good starting value 
of (S z ), one can get the renormalized locator 7 and 
the magnon spectral function D{lu') from the CPA self- 
consistent solution. Furthermore, by using the Callen 
equation (|32p . one can calculate new values of magneti- 
zation (S z ), which are re-inserted in 7. 



E. The full self-consistent loop 

In this subsection, we describe the full self-consistent 
loop for treating the FKLM in order to incorporate the 
theoretical methods mentioned in the above subsections, 
where a conducting electronic problem and a local mag- 
netic problem are solved separately. 

Such a procedure is illustrated in Fig. [T] and pro- 
ceeds as follows. For a given temperature (e.g. T = 0) 
and an initial magnetization {S z }, the electronic part of 
the model Hamiltonian @ is solved by using a retarded 
single-particle Green function within ISA method (in the 
Sec. Ill A|) . The result provides the approximate self- 
energy of the conducting electrons in the concentrated 
KLM (CKLM). In the diluted ferromagnetic systems, 











No 








+ 
H 


/ 


\ 




\ V 



IS SOLUTION OF <S Z > ) 
I SELF-CONSISTENT? j 



ELECTRONIC SELF-ENERGY 
MANY-BODY PROBLEM IN CKLM 



ELECTRONIC SELF-ENERGY 
SELF-CONSISTENT CPA IN DKLM 



ELECTRONIC PART 



MAGNON GREEN FUNCTION 
SELF-CONSISTENT CPA IN HM 
CALLEN EQUATION: <S Z > 



MRKKY 

EFFECTIVE EXCHANGE INTEGRALS 



LOCAL MAGNETIC PART 



FIG. 1: Flow chart describing the self-consistent calculation 
of the diluted ferromagnetic Kondo-lattice model. 



for a given occupation of the conduction band, the self- 
energy of the conduction electrons can be determined 
self-consistently using the coherent potential approxima- 
tion (CPA) within a dynamical alloy analogy (in the Sec. 
Ill Bp . To this step, the treatments of the electronic part 
in the diluted KLM (DKLM) is finished and one gets the 
single-electron Green function of the diluted ferromag- 
netic systems. In order to study the magnetic proper- 
ties, we use a modified RKKY (MRKKY) theory that 
results from mapping of the s — f interaction |T]) onto 
an effective Heisenberg model (HM) just as it is done 
in the conventional RKKY theory (in the Sec. Ill C[) . 
The MRKKY theory 5 ^, takes into account higher order 
terms of the induced conduction electron spin polariza- 
tion and results in an effective exchange integrals J*J 
mainly determined by the electronic self-energy. The 
latter incorporates already the dilution and disorder of 
the local magnets. After using MRKKY theory, we get 
the effective Heisenberg-exchange integrals between any 
two lattice-sites. Finally, considering random occupa- 
tion of the magnetic ions, the effective Heisenberg model 
is solved in a CPA framework self-consistently and the 
new magnetization (S z ) can be derived from the retarded 
magnon Green function and the well-known Callen equa- 
tion (in Sec. Ill D[) . The entire procedure of solution is 
repeated until the solution of (S z ) is self-consistent. Fur- 
thermore, one can change the temperature and repeat 
the entire (S z ) self-consistent calculation until one gets 



III. MAGNETIC AND ELECTRONIC 
PROPERTIES 

We have evaluated our theory for a local-moment sys- 
tem (S — 5/2) on a sc lattice with the lattice constant 
a = 1 where the width of the free Bloch-band is chosen 
to be W = 1 eV. In addition, the concentrations of the 
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nonmagnetic and magnetic ions is set to ca — 1 — x and 
cb = x. We are interested in how dilution and disorder 
of the localized moments influence the (interconnected) 
electronic and magnetic properties of the system. We 
start with the inspection of the electronic part in terms 
of the quasiparticle density of states (QDOS). 
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FIG. 2: Quasiparticle density of states p a of the conduction 
band as a function of energy in the 'diluted' Kondo-lattice 
model, in (a) and (a') for ferromagnetic saturation (T = 0) 
and two different exchange couplings J, in (b) for J = 0.7 eV 
and different temperatures. Parameters: sc-lattice, W = 1 
eV and S = 5/2. 

Fig. [3 (a) and (a') show the quasiparticle density of 
states of the conduction band at T = for two differ- 
ent exchange couplings J and for several concentrations 
x of the magnetic ions. For not too small exchange cou- 
plings J the spectrum is split into several quasiparticle 
subbands each of them with clear physical meaning. The 
low-energy as well as the high-energy bands in Fig. [2] 
(a) are due to the sites occupied by magnetic moments. 
We denote these bands in the following as 'correlated' 
bands. The middle structure is due to the non-magnetic 
sites ('uncorrelated' bands). Consequently the spectral 
weight of the uncorrelated band increases with decreas- 
ing x at the expense of the weight of the correlated bands. 

Let us first discuss the 'concentrated' case, i. e. x = 1. 
In the saturated ferromagnetic phase (T = 0), a spin- 
up (f ) electron has no chance to exchange its spin with 
the localized moments since all spins are parallel to each 
other and parallel to the electron spin aligned. The spin- 



flip terms in the interaction part of Eq. ([4]) do not work, 
only the Ising-like part (first part of the Hamiltonian 
H s f) takes care for a rigid shift by —\JS of the f QDOS 
compared to the 'free' Bloch DOS. 

The spin-down J, spectrum is more complex since a J, 
electron has even at T — two possibilities to exchange 
its spin with the ferromagnetically saturated local mo- 
ment system. It can emit a magnon therewith reversing 
its spin. This leads to 'scattering states' which built up 
the low-energy subband (~ —^JS). Because of the spin 
reversal, after which the charge carrier hops over lattice 
sites with exclusively parallel aligned local-moment spins, 
this subband occupies exactly the same energy region as 
the t subband. 

The I electron has also the possibility to exchange its 
spin by polarizing its local spin neighborhood what can 
be understood as a repeated magnon emission and re- 
absorption. This creates a typical quasiparticle, which 
is called 'magnetic polaron' and can be considered as 
a propagating electron dressed by a virtual magnon 
cloud. In the limit of a single electron in an otherwise 
empty band the magnetic polaron even gets an infinite 
lifetime^. The high-energy J. subband (~ +± J(S + 1)) 
is built up by such polaron states. Magnon emission by 
the I electron is equivalent to magnon absorption by a t 
electron, however, only if there are magnons in the sys- 
tem. In ferromagnetic saturation magnons do not exist. 
This is the reason why there is no high-energy part of the 
| spectrum at T = 0. It appears for finite temperatures 
Fig. S](b). 

The physical interpretation of the correlated quasipar- 
ticle subbands does not change for the diluted case x < 1 . 
For not too weak couplings J as J = 0.7 e V in Fig. [2] (a) 
one recognizes the appearance of an additional middle 
structure (~ 0) the spectral weight of which scales with 
1 — x. There is only a slight induced exchange splitting 
of the uncorrelated bands due to hybridization with the 
correlated subbands. For weak couplings as J = 0.1 eV 
in Fig.[^a') correlated and uncorrelated bands are mixed 
preventing a clear interpretation of the various influences. 
With decreasing x, the | and j spectrum become more 
and more symmetric since the uncorrelated bands are al- 
most unpolarized and the correlated part is merged for 
low x into the dominating uncorrelated bands. This phe- 
nomenon shows that for weak J and a given band occu- 
pation the dilution will decrease the charge carrier polar- 
ization. 

Fig. [2] (b) shows the temperature-dependent QDOS of 
the conduction band. At finite temperature, there are 
magnons to be absorbed by the f electrons and therefore 
scattering and polaron states also appear in the high- 
energy | spectrum. Consequently, the high-energy J, elec- 
tron (magnetic polaron) has now a finite spin-flip prob- 
ability to become a | electron. This spin-flip is not pos- 
sible for the saturated ferromagnetic phase since there is 
no high-energy 'scattering' part of the f spectrum. In 
the low-energy subband the electron hops mainly over 
lattice sites, where it orients its spin parallel to the lo- 
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cal moment, either without or with preceding spin-flip 
by magnon emission (absorption). With increasing tem- 
perature spin disorder of local moments grows up. The 
average distance for the propagating f electron to find 
a parallel local spin therefore becomes larger which re- 
sults in a shrinking width of the lower | subband with 
increasing temperature. Finally, the spectrum of f and j. 
electrons become symmetric at T > Tq, where spin bands 
are equally populated and spin-polarization disappears. 
The uncorrelated bands do hardly show up any temper- 
ature dependence, at least as long as they are separated 
from the correlated bands. 

Note that the KLM does not consider a direct exchange 
between the localized moments. The coupling of the local 
spins is therefore of indirect nature being mainly medi- 
ated by the itinerant charge carriers via interband ex- 
change to the localized moments. So, the possibility of 
a collective moment ordering is strongly influenced by 
the itinerant electron concentration n and by the con- 
centration x of the randomly distributed local moments. 
In case of more or less separated correlated and uncorre- 
lated subbands the effective filling of the correlated bands 
is given by n/x which turns out to be the decisive param- 
eter. Fig. [3] shows the influence of x, n together with the 
interband exchange coupling J on the effective nearest- 
neighbor Heisenberg-exchange integral at T = 0. 
The carrier-induced exchange interaction is long-ranged. 



T=0 T=0 




t i i i i i =i i i i i , i i 

0.5 1 1.5 0.5 1 1.5 

J (eV) J (eV) 



FIG. 3: Effective nearest-neighbor exchange integrals at 
T = as a function of the interband exchange coupling J and 
for various moment concentrations x: (a) carrier concentra- 
tion n = 0.1; (bl)-(b5): average number of electron per mag- 
netic ion: n/x = 0.1,0.3,0.5,0.7,0.9. Parameters: sc-lattice, 
W = 1 eV, S = 5/2. Note the different energy scales. 



However, by far the most important contribution stems 
from the nearest-neighbor Heisenberg-exchange integral 
Ji . It is obvious that the interplay between J, n and 
x does lead to a rather complicated behavior of the ef- 
fective RKKY exchange integrals far beyond the expec- 
tation based on conventional RKKY. 

In Fig. [3J (a), we investigate the J dependence of the 
effective Heisenberg exchange for a given carrier density 
(n = 0.1). For small interband exchange coupling J — > 
the itinerant charge carriers move almost freely in the en- 
vironment of local moments. Consequently, J\ shows a 
conventional RKKY behavior oc J 2 as can be found 
by perturbation theory. For small J and strong dilution 
(small x) the correlated band states are deeply within the 
uncorrelated band, very close to the band centre. There- 
fore, for small carrier concentrations only uncorrelated 
states are occupied giving rise to 'conventional coupling 
behavior'. The J region of 'conventional coupling behav- 
ior' is the larger the stronger the moment dilution, be- 
cause for lower x the system is more 'free-electron- like', 
i.e. less polarized (see Fig. [5] (a)). This explains why 
the observed maximum increases with decreasing 

x. For increasing J, the coupling H s f eventually starts 
to influence decisively the transport of itinerant carriers, 
e.g. by polarizing the spin of the conduction electron. 
Conventional RKKY is no longer valid. 

With J — > oo, the effective exchange integrals Jf^ sat- 
urate. Strong coupling J leads to a large gap between the 
separated correlated and uncorrelated bands (see Fig. ^ 
and a further increase of J will not change remarkably 
the physics of the KLM. For small enough particle con- 
centrations the electron hopping between two magnetic 
sites now happens no longer via the uncorrelated bands 
as in the weak-coupling case, but almost exclusively via 
the correlated band. On the other hand, the latter are 
getting narrower with increasing dilution, i.e. decreas- 
ing x. The electron mobility decreases and therewith the 
carrier-induced coupling between the localized moments. 
This may explain why the saturation value of Jf * is 
getting smaller with increasing dilution. 

In the intermediate J region the problem becomes 
rather sophisticated because of the competition between 
the different coupling mechanisms working for, respec- 
tively, weak and strong interband exchange J. 

There is no difference Figs.|3j(a) for n = 0.1 and Figs. [3J 
(bl) for fixed n/x — 0.1 for the concentrated case (x = 1). 
However, dilution [x < 1) in Figs. [3] (bl) means that 
there are less itinerant electrons to provide the coupling 
between the magnetic moments. Different from the situa- 
tion in Figs.[3j(a), explained above, the effective nearest- 
neighbor coupling decreases in the weak coupling 
region (J small, conventional RKKY) with increasing di- 
lution. The strong-coupling behavior is to be understood 
as that in Fig. [3J (a). 

Figs.[3j(b2)-(b5) demonstrate the change of with 
increasing n/x. In the strong coupling region, n/x is 
just the effective occupation of the correlated low-energy 
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subband. The disordered KLM is in its double exchange 
limit favoring parallel spin alignment because that in- 
creases the electron hopping probability. Thus, the sys- 
tem always has ferromagnetic J-f with the highest sat- 
urated value at n/x ~ 0.5. The correlated f subband is 
just half-filled, a very convenient situation for ferromag- 
netism in the KLM—. The decrease of the strength of 
Ji with decreasing x for fixed n/x is again due to the 
accompanying decrease of the free charge carrier concen- 
tration. The maximum at n/x ~ 0.5 in the strong cou- 
pling J-f -curves is the result of two competing mech- 
anisms. At first, increasing n/x for fixed x means that 
the system gets more itinerant charge carriers, so that the 
strength of tends to increase. On the other hand, 
increasing n/x means also that the correlated subband 
fills up, so that the quasiparticle hopping is going to be 
blocked due to Pauli principle. Accordingly, the tendency 
to antiferromagnetic moment ordering and respective an- 
tiferromagnetic correlations come into play suppressing 

jeff 

The competition between the two counteracting ten- 
dencies become even more evident in the intermediate 
coupling range (J = 0.25---0.75 eV). The correlated sub- 
bands are already split off especially for high concentra- 
tion x. They arc effectively half-filled, i.e. the low-energy 
part is occupied, the high-energy part is empty, when n/x 
approaches unity. We know from the 'normal' KLM that 
the antiferromagnetic phase as well as phase separation 
become likel y 54 ' 55 due to the Pauli principle. This man- 
ifests itself in negative exchange couplings between the 
localized moments (Figs. [3] (b4),(b5)). However, antifer- 
romagnetism is very sensitive to local-moment dilution 56 . 
In the diluted case (x < 1), due to the hybrid of the 
correlated band and the uncorrelated band, the itiner- 
ant carriers can enter the uncorrelated band although 
n/x = 1, i.e. they are no longer forced to hop exclusively 
over magnetic sites. That weakens the Pauli blocking 
and enhances mobility of the spin polarized charge carri- 
ers. That is the reason why J% changes from negative 
to positive values when x decreases from 1.0 to 0.6. 

Figs. 2] (Al),(Bl) show the distance dependence of the 
effective exchange integral J e ^ . The magnetic neigh- 
bors of a given magnetic ion are considered as ordered 
in 'shells'. The larger the shell number L the larger the 
distance to the given magnetic ion. The L-th shell is 
built up by L-th next-nearest neighbors. Since the effec- 
tive exchange integrals are strongly carrier-concentration 
dependent the exchange integrals are calculated for such 
n- value for which the system has the maximum Tq, i.e. 
n = n(T = Tff ax ) (see Fig.©. The inserts (a) and (b) 
in Fig. [4] demonstrate that the effective nearest- neighbor 
is much larger than (L > 2). Comparing 
part (Al) and (Bl) of Fig. 2] one recognizes that J e ^ 
originating from the weak-coupling interband-cxchangc 
(small J) is oscillating and long-range corresponding to 
the well known conventional RKKY-behavior. For strong 
J the effective exchange J e ^ turns out to be substan- 



J=0.9 eV J=0.5 eV 




T (K) T (K) 

FIG. 4: (Al), (Bl): Effective exchange integral J eff at T = 
as a function of No. of shell (L) for various concentrations 
x. The L-th shell is built up by L-th next-nearest neigh- 
bors. (A2), (B2): Temperature-dependence of the effective 
nearest-neighbor exchange integral J^" for various moment 
concentration x and two different values of J (0.9 eV and 
0.5 eV). (A3), (B3): Magnetization as a function of temper- 
ature T for the same concentration x and the same J in the 
other parts. Parameters: sc-lattice, W = 1 eV, S = 5/2, 

„ (rr\ rrimax\ 

n = n(T = T c ). 



tially stronger damped. The magnetism of the FKLM 
is then sufficiently well described by an effective Heisen- 
berg model with short or intermediate range exchange 
integrals J e ^ taking into account the first several shells. 
For decreasing x, i.e. increasing dilution, J e ^ has less 
oscillating behavior. For strong dilution x = 0.1, as 
realistic for real diluted ferromagnetic semiconductors 
(Gai~ x Mn x As), all interactions are essentially ferro- 
magnetic over all distances and are also exponentially 
damped with respect to distance. This agrees with re- 
cent ab-initio calculations 4 ^. 

Fig. H] (A2), (B2) show the temperature dependence 
of the effective nearest-neighbor exchange coupling . 
The T dependence is not very striking except for a slight 
enhancement near T — Tq that leads to more or less 
linear behaviors of the magnetization (S z ) for T — > Tc 
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(Fig. [4] (A3), (B3)). However, we do not find the remark- 
ably anomalous magnetization behavior as function of 
temperature as predicted by Sarma et. a/.— and also by 
reference 6 . Furthermore, dilution suppresses J{ ff for all 
temperatures, and its influence obviously increases with 
increasing interband exchange J. 

Fig. [5] shows the Curie temperature Tc as a function 
of the average number of electrons per magnetic ion n/x 
for various exchange couplings J and different magnetic 
ion concentration x. In Fig. [5] (a), some general features 
are worth to be mentioned: (1) Extremely low carrier 
concentrations are already sufficient to induce ferromag- 
netism. (2) The transition temperature exhibits a strik- 
ingly nonmonotonic behavior with a distinct maximum. 
This points to the importance of compensation effects in 
real systems for getting optimal Curie temperatures. The 
number of free carriers should be smaller than the num- 
ber of magnetic ions at least for not too low x. (3) There 
is always a break-down of the Tc curves at a critical 
n/x. For higher values of n/x the system gets negative 
magnon energies indicating an instability of the ferro- 
magnetic state against antiferromagnetism or spin glass 
behavior, especially for the weak J. (4) For a given n/x 
and a given J the Curie temperature decreases with in- 
creasing dilution, whilst the ferromagnetic n/ir-region in- 
creases as is also found by K. Aryanpour et. al. 30 . Note 
that the effects due to dilution and disorder, respectively, 



(a) (b) 




n/x n/x 



FIG. 5: Curie temperature Tc as a function of the ratio of 
the band occupation n and the magnetic ions concentration x 
for various interband exchange couplings J and magnetic ions 
concentration x. Parameters: sc-lattice, W = 1 eV, S = 5/2. 



are strongly interconnected and are not separated in this 
paper. In our previous paper—, we have demonstrated 
with a different method that dilution appears to be more 
important for the ferromagnetism than disorder. 

Part (b) of Fig. [5] shows Tc as a function of n/x for 
various x and J = 0.4 • • • 0.9 eV. The general features 
are the same as in the left part. However, it is inter- 
esting to observe that for the non-diluted case x = 1.0 
larger exchange coupling J leads to higher T c , as intu- 
itively expected. For a moment concentration of about 
x = 0.5 the Tc-curves appear to be almost indepen- 
dent of J, while for strong dilution x — 0.1 stronger 
J lead to lower transition temperatures. This is a di- 
rect consequence of the nonmonotonic behavior of the 
effective exchange coupling between the randomly dis- 
tributed localized magnetic moments, discussed above. 
In the highly diluted case x = 0.1 and for strong cou- 
plings J the low-energy correlated band is well separated 
from the broad uncorrelated band being itself rather nar- 
row. The 'free' carriers are therefore rather immobile 
since the chemical potential lies within the low-energy 
quasiparticle subband. The carrier immobility leads to 
a short-range J e ^ in the strong J case. For the strong 
dilution x ~ 0.1, however, the average distance between 
the local moments is large and therefore the short-range 
J e ff has only small influence on Tc- This is different for 
smaller J because the correlated band then melts into the 
uncorrelated one. That means a distinctly higher mobil- 
ity of the charge carriers which mediate the indirect cou- 
pling between the local moments and therefore leads to a 
long-range J e ^ . The long-range J e ^ will effectively en- 
hance Tc- So we understand that counterintuitively the 
lower J gives higher Tc than the larger J in the strongly 
diluted limit. The absolute values of Tc are of course 
strongly reduced with decreasing moment concentration. 
The situation is completely different for the concentrated 
case x = 1.0 because then no uncorrelated band exists, 
where 'normal' behavior is observed, namely higher Tc 
for stronger interband exchange J. For x — 0.5 both 
trends are obviously cancelling each other. The DMFT 
calculation by Chattopadhyay et. al^ shows also a non- 
monotonic behavior of the transition temperature. In 
their calculation, the weak J broadens the ferromagnetic 
regions even up to n/x ~ 3, probably due to an artefact 
of the method. On the other hand, the ferromagnetic 
regions for the strong J are limited to n/x < 1, rather 
similar to our results. Furthermore, the observed Tc~ 
reduction for strong J in the very diluted case applies for 
n/x ~ 1 in their calculation as well. 



IV. SUMMARY 

We have performed a selfconsistent model calculation 
of the electronic and magnetic properties of diluted local- 
moment systems described by an extended (ferromag- 
netic) Kondo-lattice model (s-f model) . It was the main 
goal to treat the exchange interaction and the moment 
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disorder on the same level. The electronic selfenergy 
was derived by a Green's function formalism ( 'Interpolat- 
ing Selfenergy Approach ' ISA) previously developed and 
tested for the 'normal' KLM. The ISA had to be gen- 
eralized to the diluted case, i.e. to the situation of ran- 
domly distributed magnetic moments. For this purpose 
we applied a CPA treatment to a dynamical alloy anal- 
ogy of the diluted KLM. The resulting self-energy turns 
out to be dependent on the magnetization of the diluted 
and disordered local-moment system. This quantity was 
determined by mapping the interband exchange interac- 
tion of the KLM on an effective Hciscnbcrg model. The 
effective Heisenberg-exchange integrals appear as func- 
tionals of the electronic selfenergy getting therewith a 
distinct temperature- and carrier concentration depen- 
dence. In the last step we solved the problem of the 
diluted Hciscnbcrg spin system by a generalized locator- 
CPA. Finally we arrived at a closed system of equations 
which could be solved selfconsistently for the electronic 
and magnetic properties of the diluted local-moment sys- 
tem. The latter could be discussed in dependence of the 
moment concentration x, the concentration of itinerant 
charge carriers n, the interband exchange coupling J, and 
the temperature T. 

The electronic quasiparticle energy spectrum is 
strongly influenced by the dilution and disorder of the 
local-moment system. For strong enough exchange J and 
finite temperature T the quasiparticle density of states 
consists for each spin direction of three subbands. Two 
of them ( 'correlated subbands ') are due to the existence of 
magnetic moments and are responsible for characteristic 
correlation effects which manifest themselves, e.g., in a 
remarkable temperature dependence. The third ( 'uncor- 
rected') subband, which exhibits hardly any observable 
T-dependence, is caused by the nonmagnetic sites. For 
weak couplings the three subbands are hybridizing. Since 
quasiparticle density of states and electronic selfenergy, 
respectively, decisively influence the effective Heisenberg 
exchange integrals, the special electronic quasiparticle 
structure determines the magnetic properties, too. 

The effective Heisenberg exchange integrals reproduce 
in the weak coupling limit (J — > 0) and for not too strong 
dilutions the well known RKKY-behavior (J^ cx J 2 ), i.e. 
long range and oscillating. This changes drastically, how- 
ever, for moderate and strong couplings, where due to 
electronic correlation effects a strongly non-monotonic 
J-dependence does appear. When the average num- 
ber of electrons per magnetic ion n/x approaches 1, the 
dominating effective nearest-neighbor exchange integral 
J?ff assumes negative values indicating antiferromag- 



netic tendencies. It is worth to be noted that the ap- 
proach covers RKKY as well as double exchange behav- 
ior. 

The Curie temperature Tc of the diluted spin system 
has been calculated as a function of the average num- 
ber of electrons per magnetic ion n/x for various ex- 
change couplings J and different magnetic ion concen- 
trations x. We found that extremely low carrier concen- 
trations are already sufficient to induce ferromagnetism, 
where the transition temperature exhibits a strikingly 
non-monotonic behavior with a distinct maximum. This 
demonstrates the importance of compensation effects in 
real systems for getting optimal Curie temperatures. The 
number of free carriers should be smaller than the num- 
ber of magnetic ions at least for not too low x. The re- 
sults also show that there is always a break-down of the 
Tc curves at a critical n/x, indicating instabilities of the 
ferromagnetic state against antiferromagnetism or spin 
glass phases. For a given n/x and a given J the Curie 
temperature decreases with increasing dilution, although 
the ferromagnetic n/x-region simultaneously broadens. 
As far as the dependence of Tc on J (for a given n/x) is 
concerned we observe two regimes. For high x a higher 
J means an increased Curie temperature. However, for 
lower x an increasing J suppresses Tc- 

Our investigation of the mutual influence of magnetic 
correlations and disorder effects in diluted local-moment 
systems clearly indicates for all interesting parameter 
constellations the necessity to treat both phenomena on 
comparable theoretical level, and that beyond conven- 
tional perturbation or mean field theory. The complex- 
ity of the problem did not yet allow us to include the 
inspection of antiferromagnetic phases which is intended 
for the near future. To tackle real diluted ferromagnetic 
semiconductors such as Ga\- x Mn x As it is also planned 
to combine our model study for realistic parameters with 
an 'ab initio' band structure calculation to allow for a 
direct comparison with experimental data. 
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